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This  report  is  Part  I of  the  Final  Technical  Report  of  a study 
being  conducted  by  the  Electrical  Engineering  Department  under  the 
auspices  of  the  Engineering  Experiment  Station  of  Auburn  University. 
This  technical  report  is  submitted  toward  fulfillment  of  the  require- 
ments prescribed  in  AFOSR  Contract  F44620-76-C-0054. 
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Transient  Analysis  of  a Finite  Length  Cylindrical  Scatterer 
Very  Near  a Perfectly  Conducting  Ground 


T.  H.  Shumpert 
Auburn  University 
Auburn,  Alabama 


ABSTRACT 

In  attempting  to  model  and  predict  the  magnitude  of  the  surface 
currents  induced  on  aircraft  in  the  ground-alert  mode,  it  is  necessary 
to  examine  the  effects  of  the  near  proximity  of  the  earth's  surface. 

For  thin  cylindrical  scatterers  sufficiently  far  removed  (several  wave- 
lengths) from  the  surface,  these  effects  may  be  taken  into  account  with 
filamentary  currents  on  the  scatterer  and  its  image.  However,  if  the 
scatterer  is  moved  very  near  (a  fraction  of  a wavelength)  to  the  ground, 
the  assumption  of  filamentary  currents  is  invalidated.  In  this  note  a 
transmission  line  mode  approximation  is  used  to  model  the  circumferential 
variations  of  the  surface  current  induced  on  a finite  length  cylindrical 
scatterer  very  near  a perfect  ground.  This  solution  is  compared  to 
previous  solutions  based  on  filamentary  currents.  The  results  give  clear 
indications  as  to  v/hen  the  more  sophisticated  approach  should  be  used  to 
obtain  valid  solutions  to  the  scattering  problems  of  this  type. 
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I.  iirrRODUCTioii 


Previous  investigations  have  considered  the  interactions  of  thin 
cylinders  with  an  electromagnetic  pulse  over  perfectly  conducting 
grounds  [l-5]*  Limitations  imposed  by  the  so-CEilled  "thin-wire" 
assiunptions  and  approximations  are  inherent  short-comings  [6-8].  In 
general,  these  approximations  can  be  divided  into  three  areas: 

(1)  current  is  assumed  to  flow  only  in  the  direction  of  the  wire  axis, 

(2)  boundary  conditions  are  applied  only  to  the  axial  component  of  the 
electromagnetic  field  at  the  wire  surface,  (3)  current  and  charge  den- 
sities are  approximated  by  filaments  of  current  and  charge  on  the  wire 
axis  [9-11].  The  emphasis  of  this  investigation  is  on  the  last  of 
these  approximations.  The  first  assumption  ignores  the  induced  current 
in  the  circumferential  direction,  which  is  an  appropriate  approach  [T]» 
[12]  provided  the  length  of  the  cylindrical  scatterer  is  much  greater 
than  its  radius.  V/ith  this  restriction,  the  scattered  field  is  deter- 
mined primarily  by  the  longitudinal  component  of  the  current  so  that 
the  significance  of  the  circiamferential  component  is  minimal.  It  is 
veil  known  that  for  infinitely  long  cylinders,  the  axial  component  of 
the  Incident  electric  field  produces  only  currents  in  the  axial  direc- 
tion and  the  component  of  the  incident  electric  field  in  the  circumfer- 
ential direction  results  in  only  circumferential  currents  [l3-2l].  For 
finite  length  cylinders,  either  component  of  the  incident  electric 
field  excites  current  in  both  the  axial  and  circumferential  directions 
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(22-251 . The  second  "thin-vire"  approximation  does  not  take  into  ac- 
count that  portion  of  the  axial  current  contributed  by  the  circumferen- 
tial component  of  the  incident  electric  field.  The  axial  current  caused 
by  the  axial  component  of  the  incident  field  is  much  more  sicnificant 
than  that  resulting  from  an  incident  electric  field  with  a circuraferen- 
tial  component.  This  restriction,  that  the  cylindrical  scattcrer  be 
thin,  makes  this  approximation  very  reasonable  [26].  Representing  the 
current  and  charge  densities  induced  on  a thin  cylindrical  scatteror 
by  filaments  of  current  and  charge  on  the  cylinder  axis  is  in  effect 
assuming  that  their  circumferential  variations  are  uniform  [l-l*].  This 
is  well  founded  for  a thin  cylindrical  scatterer  many  radii  away  from 
the  ground  plane  [6],  [27],  but  certainly  not  correct  when  the  cylin- 
drical scatterer  is  positioned  near  the  ground  plane  - on  the  order  of 
a radius  away.  In  this  analysis,  the  circumferential  behavior  of  the 
induced  currents  on  a thin  cylinder  is  taken  into  account  when  the 
scatterer  is  near  the  ground  plane. 

A Pocklington  type  integro-dif ferential  equation  [ll]  is  formulated 
for  the  current  induced  on  the  thin  cylinder  and  its  image  in  terms  of 
a complex  frequency.  This  equation  is  reduced  to  a system  of  algebraic 
matrix  equations  through  application  of  the  method  of  moments  [9],  [ll], 
(28-29].  This  transient  analysis  problem  employs  the  singularity 
expansion  method,  which  was  formalized  and  discussed  in  general  by 
Baum,  among  others  [30-33],  and  practically  demonstrated  by  several 
researchers  [l-3],  [3^J-37]»  The  complex  natural  resonances,  natural 
mode  vectors,  and  normalization  coefficients  are  calculated  and  com- 
pared to  those  found  through  enforcement  of  all  the  "thin-wire" 
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assumptions.  As  the  cylinder  approaches  the  ground  plane,  the  trajec- 
tories of  certain  singularities  are  presented  and  discussed.  Induced 
currents  are  calculated  for  various  geometries  and  incident  fields. 
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II.  THEORY 


Intcf'i'O-dif fevential  Equation 

Consider  a finite  length,  infinitely  thin-walled,  perfectly  con- 
ducting, right  circular  cylinder  as  shown  in  Figxxre  2-1.  The  cylinder 
is  near  and  parallel  to  an  infinite,  perfectly  conducting,  ground  plane. 
As  indicated  in  the  figure,  the  cylinder  is  of  length  L,  radius  a,  and 
height  h above  the  ground  plane.  A combined  cartesian  and  cylindrical 
coordinate  system  is  centered  on  the  cylinder  as  shown.  The  system, 
consisting  of  the  cylinder  and  the  grovind  plane,  is  illuminated  by  a 
transient  incident  field  of  electromagnetic  radiation.  The  incident 
field  is,  by  definition,  that  field  which  would  exist  if  the  cylinder 
and  groimd  plane  were  absent.  As  shown  in  Figure  2-2,  the  incident 
electromagnetic  field  propagates  in  a general  direction  described  by 
the  angle  0^  with  respect  to  the  z axis  and  the  angle  with  respect 
to  the  X axis.  It  is  then  desired  to  obtain  the  induced  currents  on 
the  cylinder  as  a function  of  time. 

By  application  of  image  theory  [l8],  [38-39],  the  cylinder  and 
incident  TEl  transient  plane  wave,  in  conjunction  with  the  perfect 
ground  plane,  are  transformed  into  an  equivalent  problem  consisting  of 
the  "original  cylinder"  - to  be  called  the  object  cylinder  - and  its 
image  - to  be  named  the  image  cylinder.  According  to  image  theory,  the 
incident  field  must  be  imaged  also  - producing  the  equivalent  twe-body 
pr6bl<‘m  shown  in  Figure  2-3.  Two  individual  coordinate  syster.s  are 
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IMAGE 

CYLINDER 


OBJECT 

CYLINDER 


(Note:  due  to  image  theory,  = 0^  and  = <|>^  ; subscripts  refer  to 
coordinate  dependence  and  superscripts  refer  to  incident  field 
dependence. ) 

Fieure  2-3.  Equivalent  Inace  Theory  Problem 


defined  by  Figure  2-3,  indicated  by  the  subscripts  "o"  and  "i",  repre- 
senting the  object  and  image  coordinate  systems  respectfully.  Note  the 
redundancy  of  defining  a zq  and  zi  axis  since  the  cylinders  are  iden- 
tical and  parallel.  With  respect  to  the  electromagnetic  excitation, 
the  term. "incident  field"  shall  now  be  understood  to  represent  the 
field  plus  its  reflection  from  the  ground  plane.  The  surface  currents 
induced  on  the  object  and  its  image  by  this  incident  TEM  plane  wave  are 
considered  as  equivalent  source  currents  radiating  in  free  space  (l8]. 
Thus,  the  principles  of  free  space  Green's  functions  may  be  used  to 
compute  the  scattered  field  at  an  arbitrary  field  point  in  space. 

Referring  to  Figure  2-^4,  define  appropriate  magnetic  vector  poten- 
tials for  the  object  and  its  image  as  follows: 

Aq(Rq)  = the  magnetic  vector  potential  of  the  object  in 
object  coordinates 

Ai(Ri)  = the  magnetic  vector  potential  of  the  image  in 
image  coordinates 

Rq  = a general  field  point  in  space  measured  from  the 
object  coordinate  system 

I 

Rq  = a general  source  point  on  the  object  cylinder  with 
respect  to  its  coordinate  system 

^ I 

Hi  ®=  a general  source  point  on  the  image  cylinder  with 
respect  to  its  coordinate  system 

RjL  •*  the  same  general  field  point  in  terms  of  the  image 
coordinates 

Thus,  the  magnetic  vector  potential  of  the  object  is  expressed  as 


So(Ro)  = ^ 


Go(Ro»^o  ) 


(2.1) 


where  primed  indicates  source  points,  unprimed  indicates  field  points. 
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Figure  2-h.  General  Source  and  Field  Points 


and 


Go(Ro;Ro  ) = "the  free  space  Green's  function  in  object 
coordinates 

) = the  surface  current  density  radiating  in  free 
space. 

This  Green's  function  has  the  general  form  of 


Go(Ro'»^o'^  = 

I V^o'  I 


(2.2) 


Implicit  in  this  equation  is  the  assumption  that  the  temporal  variation 
of  the  fields  is  e®^,  where 


8=0  + Ju), 


(2.3) 


the  complex  frequency  variable,  with 


Y = s/c 


c = the  speed  of  light  in  free  space. 


(2.1;) 


From  Figure  2-l»,  define  a cylindrical  coordinate  system  superimposed 
upon  the  cartesian  coordinate  system  in  the  usual  manner.  Through 
simple  geometry, 

f * * *Pl/P 

|Ro-Ro  I “ fpo  +Po  -2poPo  cos  (<j>o-'*'o  ) + (z-z  ) ] ' (2.5) 


and  for  this"  circular  cylinder,  Pg  = a. 
Therefore, 

Aq(Rq)  = Aq  (pq,  4>o  » z) 

U /*  _ til 

w-w/7 


(2.6) 


t t t 


o Go(po»4>o»5';Po  ❖o  ,z  )a  dz  (2.7) 


_ I*  _2it 

r r . . 


“ liw  / / )Gg(Po»^o.z;4'o  »z  )a  '^'{•o  <^2 

o ^ o 

(2.0) 

y 

^(Po>'f’o»z)  ~ J J ^o(*f’o'»2  )e  d(»r>  dz  , (2.9) 


• • t I 


o ' o 


Ri 


where  = [Pq2  + - 2pq  a cos  (<}'o-<r>o  ) (z-z  • (2.10) 

Upon  acceptins  the  first  "thin-wire"  approximation. 


I I 


J^(4*o  ) “ ^('t’o  ^ 


(2.11) 


such  that 


■^^(pQ  »*^0  »^  ) ~ ^^(pO>’^0»^)  ^2 


(2.12) 


and 


Ao(Po.4>o.z)  = dz^  /*  f Ko((>o',z')£;]^  ad^.o’dz',  (2.13) 

•/  o J o 


or  simply 


^Oz(Po»4’o»^)  “ ^ f f Kft(4>r>' .z*  )e  adijir.'  dz'  . 

J ° Rl 


(2.1J4) 


This  process  can  be  repeated  for  the  image  cylinder,  producing  the 
similar  equation 


V 

Ai2(Pi.«i.z)  = T;!  J J Ki ( z')£l^  ad (Ji’dz’  . 

o *0  Rg 


(2.15) 
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where 


R2  = | = [p^2+a2_2p^a  cos  )+(z-z' (2.l6) 

Up  to  this  point,  the  kernels  of  (2.lU)  and  (2.l6)  arc  exact,  in  that 
the  integrations  are  over  the  surface  of  the  cylinders.  The  current 
has  not  been  assumed  to  exist  only  on  the  cylinder  axis,  which  would 
result  in  an  approximate  kernel. 

As  pointed  out  by  others  [l-U],  [6],  the  circumferential  variation 
of  the  currents  can  be  described  as  uniform  when  the  cylinder  is  many 
radii  away  from  the  ground  plane,  but  this  approximation  becomes  poor 
when  the  cylinder  is  near  the  ground  plane.  Taylor  [liO]  has  derived 
expressions  for  the  circunferenti(il  variations  of  the  axial  current  on 
an  infinitely  long  cylinder  over  a ground  plane  in  a static  mode.  As 
pointed  out  by  Taylor,  these  resulting  equations  are  also  applicable 
to  electrically  thin  cylinders  separated  a short  distance  from  the 
ground  plane  and  to  finite  length  cylinders,  provided  the  length  is 
much  greater  than  the  height  above  the  ground  plane.  Utilizing  the 
equations  of  Taylor  [l^O], 

Ko(«o'.2')  = . (2.17) 

where 

. [1  - (a/h)2ll/2  . 

^ ~ ^ * (a/h)  cos  (2.10) 

Io(''^')  “ axial  variation  bject  surface  current. 

8e;;artibllity  of  the  current  into  it  o distinct  functional  variations 


P'  becomes  a better  approximation  as  the  cylinder  becomes  longer  with  re- 

spect to  its  radius.  Note  that  as  h becomes  leirge,  the  circumferential 
variation  of  the  axial  current  becomes  uniform,  as  desired  for  a cylin- 
der far  removed  from  a ground  plane  [l-3].  Therefore,  by  assuming  the 

r 

current  to  behave  in  this  manner,  the  magnetic  vector  potentials  are 
composed  of  exact  kernels,  in  the  sense  that  the  current  resides  upon 
the  cylinder  surface  as  opposed  to  the  cylinder  axis;  the  current  is 
uniformly  distributed  about  this  s\irface  when  the  cylinder  is  at  a far 
distance  from  the  ground  plane;  and  the  current  becomes  nonunifomly 
distributed  as  the  cylinder  draws  near  the  ground  plane.  The  results 
of  these  two  approaches  - approximate  kernel  with  loniform  circumfer- 
ential variation  of  the  axial  current  and  exact  kernel  with  an  assumed 


% 

C:,  A, 


circumferential  variation  of  the  axial  current  - will  be  examined  and 
compared. 

Returning  to  (2.17) » the  image  currents  are  similarly 


where 


or 


f.  r*  M = ri  - fa/h)2ll/2 
^ 1 + (a/h)  cos  (<^i‘+Ti)  * 


t U •)  = tl  -(a/h)2ll/2 

1 1 - (a/h)  cos 


(2.19) 


(2.20) 


(2.21) 


and 


li(z')  = axial  variation  of  image  s\irface  current. 

Note  the  difference  in  q (2.21)  and  f^  (2.l8)  due  to  the  coordinate 

references  chosen.  The  magnetic  vector  potentials  become 
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/•  L /•2it 
O / ^ 

rt  Jo  Jo  2n&  ' 


XotPo**o.z)  = fizl77 


(❖o'  ad(i)o'dz' 


OVVO 

Ri 


(2.22) 


Ai(pi 


,^i,z)  = X fo 


o •'o  2va. 


^fi((}>i'  )e~7.^£  ad(»i'dz'  . 
R2 

(2.23) 


Define  two  functions,  and  F^,  as 


such  that 


^'o^Po»'f’o»z»^')  = f f^(<i>^' )e~^^l  a d(»^'  (2.2U) 

Jo  RjL 

Fi(pi,4>i,z,z' ) = f fi  (❖■■•  ' )e~^^^  a d({n  ' , (2.25) 

Jo  R2 

J/* 

I 2^^  ^o(Po»‘i’o>^>^  (2.26) 

o 

Ai(Pi.<|.i,z)  = Aj  Fi{pi,(^i.z,z')d^'  . (2.27) 


Drawing  upon  the  principles  of  image  theory,  the  currents  on  the  object 
and  image  are  related.  At  equivalent  points  on  the  object  cylinder  and 
Its  image,  the  currents  are  equal  in  magnitude,  but  opposite  in  sign. 
Stated  simply. 


Io(*’ ) = “Ii(z' ) = I(z’ ). 


(2.28) 


Therefore 


XL  * 

tia  ^ >o(Po»4’o»^>^' )Az'  (2.29) 


T- 


(2.30) 


(p  2 ) “ “ 


l<ii  / 2va 

*^r% 


7.’ ) Fi(pi,(}>i,z,z*  )dz' . 


Locate  the  field  point  at  some  general  point  on  the  surface  of  the  ob- 
ject cylinder  as  illustrated  in  Kicure  2-5.  From  the  law  of  sines. 


Pi  sin  <J)i  = a sin  a = a sin 


(2.31) 


and  from  the  law  of  cosines. 


Pi^  = a^  + Uh^  + Ijah  cor. 


(2.32) 


Therefore,  when  Ai  is  evaluated  on  the  object  surface. 


Pi  = [a2  + >jh2  + l4ah  cos 


(2.33) 


. _ . -1  ) p.  sin  

i sin  j JjjC;  ^ 


(2.3»*) 


Note  that  when  is  evaluated  on  the  surface  of  the  object. 


Po  “ a 


(2.35) 


Thus, 


^ (Po»^o.z)  1 = Ao(<)o,z) 

JSo 


(2.36) 


^(Pi»<J>i»z)  = (<{>o»2)» 


(2.37) 


which  demonstrates  their  functional  dependence,  upon  this  evaluation. 


Therefore,  on  the  object  cylinder  sui'face 


ho 


y*^  Fo(^o»z.z')<iz'  . (2.38) 


with 


r 1 I 

Fo(4'o.z.z  \ = v/l^VlT  lT(a/h)co3  (Jq'  ^ 


(2.39) 


RjL  = [2a2_2a2  cos((Jo-<'o' ) ■*■  (z-z')^]^/^ 


(2.1^0) 


and 


with 


AiUoy^)  = -^z^ 


^ f 

•'O 


^471  2iia 


^ l(z'  ) , , / . M /I  ' 

!•  ^ ( ^>0 » z , z )d2  , 


(2.1»1) 


r 2v  -7^2 

) =y/l^/h  lZXa/h)co3  <>i ' R^  a d«i ' (2.U2) 


R^  = |2a2+l4h2+Jiah  cos  (>o‘'2[ a2+lih2+liah  cos  • 

a cos(i{>i-(fii)  + (z-z')2|l/2  (2.I43) 

j a _sin  .6.0 j 

- sin  ( [a2+l4h2+ltah  cos  <>0]^^^  ' • (2.^4l4) 


The  current,  l(z'),  induced  by  the  incident  field  produces  a scattered 
field.  The  total  field,  composed  of  the  superposition  of  incident  and 
scattered  fields,  must  obey  certain  boundary  conditions,  which  enforce 
uniqueness  [lO].  The  total  field  is  defined  as 

« incident  field  + scattered  field 

= + £«  . (2.lt5) 
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The  scattered  field  is  related  to  the  induced  currents  and  charges  by 


Es(r)  = -sAS{R)-V<>s(n), 


{2.U6) 


vhcre 


AS{R)  = the  total  "scattered"  magnetic  vector  potential 
= the  total  "scattered"  electric  scalar  potential, 
vhich  can  be  related  to  A^  through  the  Lorentz  gauge  condition. 


«®(R)  = 


7 ■ As(r) 

-WoEoS 


(2.147) 


Equation  (2.1»6)  becomes 


Es(r)  = 


= -s  A^(H)-  ^ • V[V  • AS(r)1 

UoCoS^ 


(2.I48) 


The  charge  distribution  need  not  be  knovm  since  the  charges  have  been 
related  to  the  current  through  the  Lorentz  gauge  condition.  Kote  that 
since  y = s/c  and  c = I/v/PoEq,  = s^vio^o'  Thus, 


ES(r)  = -s  [^AS(r)  - ^2  • A®(R)]  ] . 


(2.1*9) 


On  the  surface  of  the  object  cylinder,  the  boundary  condition  is 


fi  X ^ = 0, 


(2.50) 


with  fl  being  the  outward  normal  unit  vector  on  the-  cylinder  surface. 
This  boundary  condition  can  also  be  represented  by 


E^^^  ss  — 

ten  Srt  ~ tan  S 


(2.51) 


wliich  nerely  states  that  the  tangential  cemponents  of  the  incident  and 
scattered  fields  must  cancel  on  the  object  surface  in  order  to  produce 


the  appropriate  boundary  condition.  Coiiibining  (2.U9)  and  (2.51)  re- 
sults in 


Js„  = - 72  'C'  ■ Je 


(2.52) 


Since 


AS(R  ) = Ao^  + = Aoj^a-  + Aij^az  , 


(2.53) 


fo  [f’o(<>o>z»^’)-^i('^o.z.z  |] 


dz' 


(2.5!*) 


Incorporating  the  ideas  proposed  by  (2.52)  and(2.5^)»  one  arrives  at 


E 


inc 

tan 


v(v-y 


f, 


[Fo(<'o.^.^  ) - 
Fi((|.cvZ.z’)l  tiz' 

-*  (2.55) 


Since 


(-l4  7iCos)s[l-^2  =-noCoS^[l-^2 

= (V(V-)  - Y^J.  (2.56) 


(2.55)  becomes 


•'0 

[Fo(<>o>z,z' )-Fi((Jo,z.z' )]az'  . (2.57) 


llie  differential  operator  [I4I]  readily  reduces  since  it  acLs  on 


a vector  with  only  a z component. 


1 

i 
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grad  div  V = V(v*V) 


(2.5O) 


[V(V)  - y2]  = 82 


dp3z 


di^QdZ 


^ &z-y^  . 


8z 


The  boundary  conditions  implied  through  (2.57)  are  simply 


(2.59) 


j,inc 

■"z 

^inc 

'4> 


-E  on  So' 
z 


-E“  on  So 
9 


(2.60) 


Through  application  of  these  two  equations,  it  is  apparent  that  for 
finite  length  cylinders  an  nocial  component  of  current  is  created  by 
both  an  axial  and  a circuiiferential  component  of  the  incident  field 
[22-25].  Nevertheless,  accepting,  the  second  "thin-wire"  approximation, 
boundary  conditions  v;ill  only  be  enforced  on  the  axial  component  of  the 
incident  field.  This  leaves  the  integro-dif ferential  equation 

(11^  - I V)  )- 

Fi(<'o»z.z' )]^z'  . (2.61) 

Referring  to  Figure  2-6,  the  total  incident  field  can  be  formulated 
on  a general  basis.  is  shown  in  the  plane  of  the  two  descriptive 
coordinate  directions  a and  b.  An  electric  field  normal  to  this  plane 
will  not  produce  a z component.  Thus, 

2*"°  = Eio"'''  * 22  P-’"’ 


= El  exp 


+ E2  exp  I -Y 


z cos  + a cos[  (tt/2)-0^] 

t 

z cos  0^  + c cos[  (it/2)-0^] 


(2.62) 


ginc  _ I f ^ 0^  + a sin  OiJ  | 

+ E2  I 6^  + c sin  I 

= exp  I ~yz  cos  oi  - Y sin  0^  |^x  cos  <{»i  + 


y cos((J^  - Tr/2) 


■li 

+ E2  I ~y^  Oi  - Y sin  oi  1^  X cos  (ir-,^i)  + 

)i  - TT/2)jj,(2. 


63) 


y cos((>^ 

where  the  principle  of  direction  cosines  [Ji2]  has  been  fowarded. 
Simplification  gives 

gine  _ j?^g“YZ  cos  0^  -y  sin  e^[x  cos  + y sin  <J>^] 

, p -yz  cos  0^  -y’ sin  0^[-x  cos  <J)^  + y sin  (})i]  (2.6I4) 

+ bpC 

Since  only  the  axial  component  of  the  incident  field  is  to  be  used, 


^zi  ^ 1^1 1 ~ sin  Oi 

^Z2  “ “1^21  sin  0^  = -Ep  sin  0^  , 


(2.65) 

(2.66) 


such  that 


,inc 


= El  sin  O^e 


-yz  cos  0^-Y  sin  0i[x  cos  + y sin  <Jii] 


-E2  sin  0^-Y  sin  0^[-x  cos  <}i^  + y sin 

^ - (2.67) 


On  the  ground  plane,  E^  must  satisfy  the  boundary  condition 


= E*"'  = 0. 

^ (Ground  Plano  ) ^ (x  = -h ) 
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(2.68) 


t 


Enforcins  this  requirement  on  (2.67)  necessitates  that 


El  = E2e 


-2Yh  sin  6^  cos  4)^ 


(2.69) 


Defining 


Eo  = -El  , 


(2.70) 


then 


E2  = -Eoe 


2Yh  sin  cos 


(2.71) 


such  that 


,inc 


-Yz  cos  0^  cos  + y sin  ) 

2Yh  sin  0^  cos  - Y sin  oi[-x  cos  + y sin 

(2.72) 


-Eq  sin  O^e 


-c 


or 


glnc 

z = -Eq  sin  O^e 


i^-YZ  cos  0^  - Yy  sin  0^  sin  , 

r -YX  sin  0^  cos  (|>^  ^YX  sin  0^  cos  ({'I  + 2Yh  sin  0^  cos  (Jilj 
^ (2.73) 


Comparison  of  this  incident  field  to  that  of  Unashankar,  et.  al.  [3], 
is  favorable.  Letting  = l80°. 


,lnc 


1 = -Eq  sin  eie-Y^ 

= Ti  *- 

^-Y  sin  ei(x  + 2h)J  , 


(2.7>+) 


,lnc 


and  evaluating  E on  the  cylinder  axis  (approximate  kernel)  instead 
of  on  i^c  surface. 
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inc 


s 

X = 


1 -Yz  cos 
-Eq  sin  O’e 


Oi 


_ ^-Y2h  sin  elj  (2.75) 


which  is  identical  to  Umashankar' s [3]  field.  This  comparison  is  made 
note  of  since  the  induced  current  found  by  Umashankar  shall  be  compared 
to  the  current  found  through  this  exact  kernel  formulation. 

Before  evaluating  the  incident  field  on  the  cylinder  surface, 
the  integro-differential  equation  {2.6l)  needs  to  be  examined  more 
closely.  The  unknown  quantity,  l(z'),  is  not  a function  of  (Jq.  It  is 
apparent  that  by  letting  <{>o  = Oi,  any  particular  angle,  solving  for 
I(z*)  and  letting  ifo  ~ “2»  ^ different  pai'ticular  angle,  again  solving 
for  l(z')  - the  two  solutions  must  be  identical.  The  implication  is 
that  in  order  to  determine  the  unknown  induced  current,  boundary 
conditions  need  not  be  enforced  all  over  the  object  cylinder  surface, 
but  just  at  one  particular  value  of  (Jq.  The  circumferential  variation 
of  the  current  has  already  been  assumed  to  be  of  the  foim  expressed  by 
(2.IO)  and  (2.21),  such  that  a solution  for  l(z')  at  any  particular 
value  of  vill  readily  result  in  a general  solution  for  l(z',<Jo). 
Through  examination,  = 0°  seems  as  profitable  as  any  other  choiee. 
Thus,  the  integro-differential  equation  reduces  to 


(-Uireos)E 


inc 


^o 
♦o  = 00 


Fi(z,z' )ldz' 


(2.76) 
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where 


Tjzii.')  = FoUo.^.2')  ittst?; 


a d({>o ' 

^1 

ri  = Ri  I = [2a^(l  - cos  (^o')  + (z  - 

Uo  = 0° 


"I 


Fi{a,z')  = FiC^o.z.z')  = QO=  yjl-a/h 


/2it  ^ 

l-(a/h)  cos 


e~^^^  a d(}>j^' 


r2  = R2  = 0°  ~ [(2a2  + U ah)  (l-  cos  (fi')  + ^h^  h 

(z-z' )2]l/2 


Evaluation  of  the  incident  field  results  in 


E 


inc 


= E 


inc 


o 

❖o  = 0° 


y = o 
X = a 


_ . -yz  cos  0^1  -ya  sin  0^  cos  Ai 

-Eq  sin  O^e  ' le  ' ^ 

^ya  sin  0^  cos  <}i^+  2yh  sin  0^  cos  (|>ij 


By  defining 


3^  = sin  0i  cos  ({>^  , 


the  Incident  field  is  simply 


,inc 


So 

«o  = OO 


-Eo  sin  Oie-^'" 

[j_^Y2Sl(n*  h)] 


— -T-* 

cos  (^Q 
(2.7Y) 

(2.78) 

(2.79) 

(2.80) 


(2.81) 

(2.82) 

(2.83) 
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Through  oxaminaLiun  of  this  equation,  the  appearance  of  a phase 
difference  between  the  two  terms  inside  the  brackets  is  obviously  that 
caused  by  the  "original"  incident  field  and  the  same  field  reflected 
fron  the  ground  plane  arriving  at  some  later  time. 

The  integro-diffcrential  equation  (2.76)  can  be  altered  slightly 
in  form  and  notation  in  order  to  better  represent  the  problem.  Since 
the  differential  operator  does  not  operate  on  z' , 


(-Uncos)E 


inc 


z 


® J 

r iw 

-y2\  I 

So 

J 

0 2.a 

0 

0 

u 

Fi(z,z')] 

esent  the 

complex 

frequency  ( 

s)  j 

- r 

I(z'  ,s)(d^ 

J 

So 

2va 

4>o  = 

00 

Fi(z,z'  ,s 

[FqCz.z'  ) - 


(2.8U) 


Let 

and 


1 - cos  <{>0'  = 2 sin^ 


d = diameter  of  the  cylinder 

ts  2a 


(2.85) 

(2.86) 

(2.87) 


in  Fq,  such  that 


[2a2(l  - COG  (Jo')  + (z-z')2]l/2  = [d2  sin2  + (z-z')2]^/2^ 

(2.88) 

with  a similar  substitution  in  Fi. 

The  final  results  are  summarized  for  reference; 

‘L 


(-Uiieos)E^"®(s) 


«o  = 00  Jo  / 

So 

[Fq(z,z' ,s)-Fi(z,z' ,s)]  dz' 


(2.89) 
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^ If 


. = -E„(s)sin  . 


^l.e 


e-yn 


/h)cos  (fp'  r^ 


(2.90) 

a d<j)  ' (2.91 ) 


rYf-2 


F.(z.z*,s)  =Jl-(a/h)^'  —L a d^.. ' (2.92) 

^ ^ 1 * (a/h)cos  4.  ' r.  ’ 

Jq  0 c 


= [d^  sin^  1|1  + (z-z')^]l/2  (2.93) 

=([d^  + 8ah)  sin^  ^ + 4h^  + (z-z  . (2.94) 

This  intecro-differential  equation  is  to  be  solved  for  the  unknown 

induced  current  on  the  cylinder.  • 

Application  of  the  Method  of  I'or.ents 

The  integro-differential  equation  shall  be  cast  into  matrix  form 
suitable  for  a numerical  solution.  This  general  process  has  come  to 
be  known  as  the  method  of  moments  [9]>  [ll]>  [ 28-29].  Generally,  the 
"wire"  along  the  cylinder  surface  at  “ 0°  is  broken  into  segments. 
Integrals  approximated  by  the  sum  of  integrals  over  N small  segments, 
itnd  the  current  assumed  to  be  constant  over  each  individxial  segment. 

With  regard  to  Figure  2-7 » expand  the  current  in  a set  of  basis 
functions  such  that 


where 


I(z',s)  = an(s)Ijj(z' ), 
n 


(2.95) 


00(3)  « unknown  coefficient  of  constant  cui'rent  in  the 
n'th  subsection 
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SUBSECTIOIJ  EiroS 


X MATCH  POItITS 
/ 

Zm  = (m-l)A  (HATCH  POIIITS) 
ra  = 1 ,2, . . .N 


z3 


* "^2 


z”  = (n-2/3)A  (SUBSECTION  ENDS) 

n = 1,2,. . .H+1 

a'  = L/(N-1)  LENGTH  OF  ZONE 
N = NUMBER  OF  SUBSECTIONS 


z 


2 


Figure  2-7-  Moment  Method  Partitioni.ng  of  Geometry 
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In(z*)  = 


(2.96) 


11-  for  < 2.  < I 

0 elsevjliere  n = 2>3> . • • .K-lj . 

Thus,  use  is  made  of  a pulse  function  expansion  described  by  Harrinston 
[9].  This  representation  is  chosen  such  that  the  boundary  conditions, 

1(0)  - I(L)  = 0,  (2.97) 

ore  satisfied  automatically  by  allowing  the  two  end  subsections  or 
zones  to  extend  past  the  surface  of  the  cylinder  and  assujring  the 
current  on  these  zones  to  be  zero.  As  shown  on  Figure  2-T,  each  zone  is 
of  length  A,  where 

A = = length  of  a zone  (2.98) 

N = number  of  subsections  or  zones 

zn  = (n-3/2)A  n = 1,2, N+1  (2.99) 

*=  subsection  ends. 


Applying  these  concepts  to  the  integro-dif ferential  equation  (2.89) 
produces 


(-l.„Cos)  j-E„(£)  sin  ein-"'"  d (a  * 2h)j  / 

/*Z°  ^ ^ / 2 \ 

^s)_Fj^(2,z' ,d)]  dz' 


(2.100) 


Forcing  this  equation  to  be  satisfied  at  discrete  natch  points  a-mounts 
to  choosing  delta  functions  as  testing  functions  (9].  The  natch 
points  are  the  center  of  subsections  as  showm  in  Figure  2-7,  where 
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(2.101) 


Zn  *=  (in-l)A  m = 1,2,. ...N 
a natch  points. 

As  pointed  out  by  Harrington  [9]»  the  derivatives  nay  be  carried  out 
ana'’ytically  or  approximated  by  finite  difference  techniques.  Both 
avenues  of  approach  were  investigated,  with  the  decision  going  to 
finite  differences  due  to  its  ease  of  evaluation  and  simplicity.  Using 
finite  difference  approximations,  where 


il 


f ‘ 


dS  “ Az)-2F(z)  + F(z-Az))  , (2.102) 


the  integro-differential  equation  i: 


) i-Eo(s)  sin  eip-’f’'*--  h-  2h)j  | 


_ ^ °n(s ) 1 r 

^ 2na  A^ 
n 


zn  + 1 


»s)  -(y?a2  + 2)* 


?’o(zm-l»z' »s)  -Fi(zn+jL,z' ,s)  + + 2)» 

^i(^ni»2  »s)  —Fi ( , z ,s)j  dz  n — 2,3,....  H— 1 

n = 2 , 3 , • . . • >^-l  • 


(2.103) 


This  equation  is  placed  in  the  form  of  a matrix  equation. 


7(s)  = l(s)  I(s), 


(2.10!^) 


where  a single  bar  represents  a column  matrix  or  vector  and  double 
bars  indicate  a square  matrix-.  Let  the  matrices  be  defined  as 

V(s)  = the  source  vector  = [v^], 

where 
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Vm  = the  matrix  elements  of  V(g) 


(-IittCqs)  -Ko(s)  sin 


0^ r -YS^a  ^YS^(a+2h) 

L J I 


m = 2,3, . . .N-1 


(2.105) 


where 


where 


l(s)  = the  response  vector  = [in]  , 


the  matrix  elements  of  l(s) 

Of),  unknown  coefficient  of  constant  current  in  the 
nth  zone 

n = 2,3, N-1  ; (2.106) 


Z(s)  = the  impedance  matrix  = [smn]  > 


Zran  = t)ie  matrix  elements  of  Z(s) 


^5,n+l 

C^o^^m+l.z' .s)  - (y^A^  + 2)  Fo(7.n,z',s) 

“'z^i 

+ r’o(zni-l>z' >s)  - Fi(zni+l,z' ,s) 

+ (y^a2  + 2)  Fi(zm,z',s)  -Fi  ( Zj^_l , z ' , s )]  dz'. 


n = 2,3 N-1 

m = 2,3. . . .N-1 


(2.107) 


"With  simplification  in  mind,  let 


rzn+1  ^2tt 

ll0„(-m.  = ) = / J YTi^ 


e~^''l  , , 

(2.108) 


with 


ri  = (d?  sin?  V ^ (-m  - z')?]^/?  , 


(2.109) 


I 


and 


Win  Jq"'  ^ -\a/h)cor^ 


e-^^2  , , 

— d<Ji  dz  , 


(2.110) 


= [(d2  + Cah)  sin^  ^ + Uh^  + (zrn  - z’)2]l/2 


(2.111) 


Thus,  Zjnn  redefined  as 


- (y2a2  + 2)  Wj^(zi^,s)  + 

H|j^(zjn_i  ,s)  -Hi^^.( Zjjj+1  ,s ) + (y^a2  + 2)  W^(znis)  - 
W^(zm-l.s)|  (2. 


(2.112) 


The  intecrals  defined  by  (2.10S)  and  2.110)  present  problems  in 
numerical  evaluation  - but  these  can  be  overcome.  Tcsche  [^*3]  provides 
methods  of  treatins  integrals  of  this  tiT^e.  Kote  that  the  intecrand  of 
is  never  sincular  due  to  the  obvious  fact  that  all  terms  in  the 
radical  (2.11)  are  creator  than  or  equal  to  zero  except  the  term  l<h2 
which  is  never  zero.  Therefore,  numerical  evaluation  of  offers  no 
problem.  However,  Hp^(zni,s)  presents  some  problems;  in  that,  when 
4>o  = Oo  and  z'  = z^^,  the  integrand  is  singuJLar.  H£^(zni,s),  Hq^  ( Znj-i , s ) , 
and  Wg^(zim.i,E)  shall  be  exa.mincd  separately. 

At  4>o  ~ '^'We  integrand  of  Ho(zni+i,s)  is  singular  v’hen  zjn+i  = z' , 

mn 

where  z'^  ^ z'  ^ z*^"*"^.  Hence,  at  this  point. 


z"  Zra+1  1 


(2.113) 
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Utilizinc  (2.99)  and  (2.101) » reproduced  here 


Zm  = (m-l)A 

n = 1,2, N+l 

an  = (n-3/2)A 

n = 1,2, ll+l  , 

singularity  occurs  when 

(n  - 3/2)A  <_  [(m+1) 

-1]A  £ ( (n+l)  -3/2]A 

(2.111*) 

or  when 


-3/2  £ (n-n)  £ - 1/2. 

Since  m and  n are  intecerr. , "m-n"  mot  be  an  integer  between  these  two 
fractions.  Thus,  Hg^^( 2^^+^ » ^ ) is  singular  when  n-n  = -1  or  when  n = n-1. 

In  a like  manner,  Hg^(zj.,,s)  is  singular  when  <>q'  = 0®  and  Zn  ~ z'. 
Thus,  the  singularity  oecurs  when 

1 ^ 1 . (2.115) 

which  can  be  expressed  as 

(n  - 3/2)A  £ (n-l)A  _<  [(n+l)  - 3/2)A  (2.116) 


or 

- 1/2  _<  (in-n)  ^ 1/2. 

Through  the  same  reasoning,  { z^ , s ) is  singular  when  m = n. 

It  is  apparent  that  Hg^^(zjn_i  ,s ) is  singula.r  when  m = n+l. 

When  • ..use  respective  conditions  do  not  exist,  the  integrals  are  readily 
amena  to  numerical  integration  techniques,  but  at  a singular  point, 
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numerical  techniques  become  suspect  at  best.  Following  general  proce- 


dures outlined  by  Tesche  let  the  integral  of  I!£^(7.ni,s) 

singularity  be 


= K£^(ziii,s)  when  m = n , 


with  defined  by  (2.108). 


Transform  variables  in  (2.108) 


let 


♦ = <{>0 
Z = Zja  - z 
dz  = -dz' 


and  as  for  limits  of  integration 


Thus, 

/*!  /’2m  ^ ^-yr 

~ Jt^  Jo  ^ +(a/h)cos  (Jr  , 

“ 2 

with  r = [z2  + d^  sin^  ^]l/2  , 

which  is  one  of  the  singular  integrals  to  be  evaluated.  In 
manner,  let 

T2  = H£^(ziii+2^ ,s)  when  m = n-1 
ana  aft»r  transforming  by  allowing 


at  its 


(2.117) 

as  follows: 


(2.118) 


(2.119) 


(2.120) 


(2.121) 

a like 


(2.122) 
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♦ = ^o' 


dz  = -dz* 

J«  = n-l  = I , „.i  = -|  _ 

f 2 Z"^*  1 e"^’' 

’^2=/-A  / I ■;  (a/h-)cos  <,  V • f2.123) 

2 *'0 

with  r defined  as  in  (2.121).  A conparison  of  (2.120)  and  (2.123)  obvi- 
ously shows  that 


T2  = Ti. 

In  addition,  T3  could  be  defined  as  H£^(zn_2^,s)  when  m = n+1  and  tlic  end 
product  would  be  that 


Ti  = Tp  = T3. 


Therefore,  the  intecral  of  the  singular  integrand  of  Hp^(zni+i,s) 
when  n = n-1,  Hg^(zjn,s)  for  m = n,  and  ( ^la-l ♦ s ) with  m = n+1  is 


-yr 


cos  <j> 


d(j>  dz 


(2.I2I4) 


with  r given  by  (2.121). 

This  Integral  is  further  pursued  in  Appendix  A,  which  concludes  with 
the  result  that 
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In  2 - 


„ ^ -2nYA Utt 

~ [1  - (a/h)^J^/^ 

ln[h/»  + a - [(h/a)2  - l]l/2]  | 2 . 

jll  ■"  1]  • (2.125) 

The  function  is  nu.T*erically  evaluated  easily,  as  the  integrand  is  not 
singular. 

The  integro-dif ferential  equation  has  been  cast  as  a system  of 
matrices,  which  must  be  solved  for  the  unknown  induced  current. 

Application  of  the  Singularity  Expansion  Method 

The  singularity  expansion  method  (SE-I)  was  introduced  by  Bauja  [31- 
33j  and  formalized  and  applied  by  many  others  [1-3],  [30j,  [3*^-3?]  as  a 
method  for  characterizing  the  response  of  scattering  objects  v/hen  illu- 
minated by  either  transient  or  steady-state  ’ ectromagnetic  radiation. 
Before  applying  SEI  to  this  transient  prob_  a brief  survey  of  the 
general  theory  is  appropriate  - highlighting  che  areas  of  particular 
interest. 

The  complex  natural  frequencies  of  the  scattering  system,  denoted 
by  So,  are  those  such  that  (2.10h),  when  expressed  as  the  homogeneous 
equation, 

2(cc)  i (sa)  = 0,  (2.126) 

has  a nontrivial  solution  for  Ksq).  The  implication  is  that  the  deter- 
minant of  Z must  vanish  at  these  complex  natural  frequencies.  Thus, 
the  ecui.tion  for  determination  of  these  natural  resonances  of  the 


induced  current  is 


det  Z(sa)  = 0 . (2.127) 

As  required  by  well-knovm  lineax  circuit  theory  and  Laplace  transform 
theory,  these  natural  resonances  must  occur  in  the  left  half  portion  of 
the  complex  s (s  = o+Jio)  plane,  defined  in  the  usual  manner  by  Laplace 
transform  theory  [l4!4-U5].  In  addition,  these  natural  resonances  can 
not  ajipear  on  the  jw  axis  and  must  consist  of  complex  conjugate  pairs. 
Tliese  facts  are  obvious,  in  order  to  produce  real  currents  in  the  time 
domain  which  eventually  eo  to  zero  due  to  the  radiation  of  enercy.  It 
has  been  shown  [35]  that  for  bodies  of  finite  extent,  the  response  func- 
tion has  only  poles  and  no  branch' cuts.  Speculation  has  been  put  forth 
[37]  that  only  simple  poles  exist  for  perfectly  conducting  bodies.  The 
assumption  that  only  simple  poles  exist  is  accepted  witliout  proof  - 
althouch  this  has  been  substantiated  nuracricaJ.ly . Other  researchers 
have  made  wide  use  of  the  simple  pole  assumption  [1-3],  [31],  [3**-37]. 


The  matrix  equation, 
l(s)  I(s)  = V(s), 


(2.128) 


has  as  a solution 

I(s)  - Z(s)-lv(s). 


(2.129) 


Through  the  inverse  Laplace  transform,  the  current  in  the  time  domain 

is 


Kt) 


Z Ts)V(s)e®^  ds. 


(2.130) 
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vhich  could  be  evaluated  nurierically  along  the  Bromwich  path.  SE-I,  an 
extension  of  well-founded  techniques  in  classical  circuit  theory,  is 
based  upon  the  idea  that  the  tine  domain  response  of  the  induced  current 
is  determined  through  knowledge  of  the  complex  natural  frequencies  of 
the  system  in  addition  to  their  corresponding  residues.  The  time 
domain  current  is  described  by  a sura  over  all  the  residues,  times  expo- 
nentially damped  sinusoids. 

Based  upon  this  concept,  SEi  assumes  that 

Z(s)“^  = entire  functions,  etc.)  (2.131) 

S— Sq 

a 

which  is  a sujsmtion  over  all  poles  in  the  complex  s plane.  Define  Rq 
as  the  system  residue  matrix  at  the  pole  Sq.  Note  that  Ra  is  a constant 
matrix,  in  the  sense  that  it  is  not  a function  of  s.  This  residue  matrix 
is  a dyadic  and  can  bo  represented  as  the  outer  product  of  two  vectors 
independent  of  s as 

h = Ma  ^0^  . (2.132) 

t 

where  Mq  is  the  natural  mode  vector  and  is  a solution  to 

Z(so)M„  = 0 (2.133) 

and  is  the  coupling  vector,  vhich  satisfies  the  equation 
Z(na)'^  Ca  = 0 

(T  denotes  transpose).  For  the  electric  field  formulation,  where 
symmetric  matrices  are  encountered,  these  two  vectors  are  identical 


Thus, 


R = C C T 

«a  *^0  '-o 


(2.13^*) 


Progrescing  one  step  further,  let  Cq  be  normalized  such  that  its  max- 


imum element  is  real  and  equal,  to  unity. 


Co  tCa  normalized] 


and  define 


(2.135) 


Cq  normalized  = Cao> 


suck  that 


^^o  “ Coq  Coq"^  » 


(2.136) 


where  is  the  normalization  coefficient. 


Since 


= E 


(2.137) 


the  normalization  coefficient  is  determined  through  consideration  of 
particular  singularity  [2],  s©,  such  that 


5p'fz(s)Z(s)“^Cp  = Cp^UCp  = CpTCp 


(2.138) 


(Note:  0 is  the  identity  matrix.  ) 


CpTZCs)  E Ba  ^5?^  Cp  = CpTcp 


E B. 


(2.13?) 
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Noting  that 


CpTz(Ep)C<j  = Cp'^Z(sp)'rTCa  = (Z(sp)TCp)Tc^  = 0, 
C^T[|(s)  _ z(sn)1C^S^-CT 


- s - Sex  " S S • 


(2.IU0) 

(2.ll»l) 


By  definition. 


lim  Z(s)  = Z(sr)  _ X , 

B*sp  — ns:: — ~ ^ '^p'  'p“  ’ 


(2.1U2) 


where 


Z' (sp) 


= ds  s = 


and  6pa  is  the  Kronecker  delta  function.  Taking  the  limit  as  s 
approaches  Sp  in  (3.1**l)  results  in 


B^Cp^z' (sp)CpCpTCp  = CpTCp  , 


(2.in3) 


from  which  we  conclude  that 


Bn  - 


CpTz'(sp)Cp 


(2.IU) 


Therefore , 


i{t) 


= — / I V C ds 

2TTj  / . ^ 

Oo-joo  a 


■0  s-S(x 


where 


n _ = B„C  V 
a a a a 
0 


(2.145) 

(2.146) 


Is  defined  as  the  coupling  coefficient.  With  the  matrix  elements  of 
V(s)  defined  by  (2.105),  let  the  incident  field  be  a step-function 
plane  wave,  such  that 


E„(s)  . E^/s. 


The  matrix  elements  of  V(s)  are  then  defined  as 
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i<- 


(2.lli7) 


1 T>  i fti  cos  9^1  -yB^a  YB^(a+2h)l 

Vj,  *=  UnCoEo  sin  O^e  ' “ Le  ' -e ' J 

m = 2,3, . . .N-1. 

Evaluating  (2.lh^)  through  the  residue  Dieorem  v/ill  produce  appropriate 
Heaviside  functions,  vhich  are  viewed  as  enforcing  causality  [37]»  [1-3). 
The  exponential  dependence  of  (2. 1^*5)  is  expressed  by 

v„0St  = De-l'‘==i  '='>=  !>‘[^-Y3ia_^YBl(a  . a.)J  (s.ne) 

where 


D = UttCoEq  sin  0^. 


(2.119) 


(2.1*48)  can  be  written  as 


vmc 


St 


= d[c'^  - 


cos  pi  + 

^s[t  - cos  pi  - P^(n  + 2h)  jJ  (2. 


150) 


such  that  (2.IU5)  becomes 


rOo+J~ 

I'*)  J Z V„(s)  da 


(2.151) 


with  Va(s)  a vector  with  matrix  elements  defined  by  (2.I50).  In  other 
words , 


Va(s)  = [vj^cst]  = v(s)est. 

For  the  sake  of  simplicity,  let 
_ _ Zm  cos  pi  + B^a 


(2.152) 


(2.153) 
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and 


T2 


Zm  COS  - 6i(a  + 2h) 
c 


(2.15‘t) 


such  that  the  matrix  elements  of  Va(s)  are 

v^e^t  = D [e^^^-^lLe^^^-^2)]  (2.155) 

Since  the  intccrand  of  (2.151 ) has  only  poles  at  s = s^,  evaluation 
throuch  the  residue  theorem  produces 


where 


i(t)  - (t ) 

a 


(2.156) 


Va(t)  = [vmci(t)] 


(2.157) 


with 

m=2,3,...N-l  (2.158) 

Note  that  the  complex  natural  frequencies,  natural  mode  vectors, 
and  normalization  coefficients  are  not  functions  of  the  incident  field. 
Only  Va(t)  is  altered  upon  a chance  in  the  angles  of  incidence.  There- 
fore, once  Sq,  3a » ^.nd  5^^  are  found  for  a particular  L,  h,  and  a, 
the  current  excited  by  any  incident  field  is  easily  found.  This  is 
perhaps  the  greatest  utility  of  SEM. 

The  singulai'ity  expansion  method,  applied  to  this  transient  elec- 
tromagnetic problem,  produces  a system  of  matrices,  which  are  solved  for 
the  i.idiccd  current  on  the  object  cylinder  as  a function  of  time. 
Heuristic  study  into  the  nature  of  those  complex  natural  frequencies. 
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natxiral  mode  vectors,  and  normeilization  coefficients  is  justified 
tlirough  a better  understanding  of  the  two-hody  transient  scattering 
Ijroblem. 

Approximations  and  Linitations  Imposed 

The  approximations  and  assumptions  which  are  accepted  throughout 
this  tlicsis  are  of  prime  importance.  VThethcr  analytical  or  numerical 
solutiojis  to  problems  are  exercised,  final  results  must  be  evaluated  in 
light  of  limitations  placed  upon  their  very  existence.  Therefore,  this 
section  brings  together  all  of  these  approximations  and  assumptions 
with  the  purpose  of  examination  and  evaluation. 

The  assumptions  and  approximations  are  as  follows: 

a)  Current  is  assumed  to  flow  only  in  tlie  direction  of  the 
cylinder  axis  (2.11). 

b)  Boundary  conditions  are  enforced  only  to  the  axial 
component  of  the  tangential  electric  field  (2.6l). 

c)  End  caps  on  the  cylinders  arc  ignored  (2.8). 

d)  The  current  is  assumed  to  be  separable  (2.17). 

e)  Taylor's  [i*0]  equations  are  adapted  for  the  circumferen- 
tial variation  of  the  currents  *(2.l8). 

f)  The  moment  method  introduces  an  approximate  numerical 
solution  (2.95). 

g) ^  The  exponential  e is  expanded  in  a two  term  Taylor 

series  for  the  solution  of  T (AI.5). 

As  previously  outlined  in  the  introduction,  the  first  two  assumptions 
restrict  the  cylinder  to  be  thin,  L >>  a.  Ignoring  the  endcaps  amounts 
to  letting  the  cross  sectional  surface  area  be  small,  a <<  X,  so  that 
ally  current  induced  on  the  endcaps  will  net  significantly  contribute  tc 
the  scattered  field  [^6).  By  assuming  t!:e  current  to  be  separable,  the 
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restriction  that  L >>  a is  again  necessary,  since  this  approximation  is 
poor  only  near  the  ends  of  the  cylinders. 

Up  to  this  point,  the  most  restrictive  clause  placed  upon  this 
formulation  is  that  the  length  of  the  cylinder  must  he  much  greater 
than  its  radius.  Thus,  the  cylinders  must  he  thin. 

Taylor  [l40]  states  that  the  equation  for  circumferential  variations 
of  the  axial,  cui'rent  (2.18)  is  applicable  to  electrica,lly  thin  cylind- 
ers separated  a short  distance  from  the  ground  plane  and  provided  the 
length  is  much  greater  than  the  height  above  the  ground  plane,  the 
expressions  can  be  used  for  finite  length  cylinders.  Concerning  the 
separation  from  the  ground  plane,  this  restriction  is  linked  to  approx- 
imating the  two  parallel  cylinders  as  supporting  a transmission  line 
m-.ii.e  (TElJ).  Typically,  this  TE.M  mode  requires  that  h <<  X,  a <<  X, 
and  >>  h.  Nevertheless,  note  that  when  a <<  h,  (2.18)  approaches 
vinity  and  the  circumferential  variations  of  the  axial  current  becomes 
uniform  - as  would  be  the  case  for  a thin  cylinder  far  from  the  ground 
plane  [l-M*  The  current  is  uniformly  distributed  in  (}>  at  a far  dis- 
tance from  the  ground  plane  and  only  becomes  nonuni fomly  distributed 
as  the  cylinder  draws  near  the  ground  plane.  Therefore,  vdien  (2.l8)  is 
actually  affecting  the  distribution  of  the  currents,  these  typical  TEl 
mode  restrictions  are  satisfied.  More  importantly,  when  the  cylinder 
is  far  removed,  circumferential  variation  of  the  currents  is  uniform, 
(2.18)  doesn't  significantly  affect  the  equations,  and  the  typical 
transmission  line  restrictions  need  not  be  satisfied. 

With  reference  to  (AI.5),  the  Taylor  series  expansion  of 
about  r = 0 is  truncated  after  two  terms: 
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- 1 _ yr  + + . . . (-1  ^ 

2!  3!  ■ ■ 


n! 


^ n! 


n-o 


where 


= (z^  + d2  r>in2  4>/2]^/^  -—  < z < ^ 


The  maximum  of  yr  is 


max 


of  yr  = y yj?-^-  + u^'  , 


where 


d = 2a. 


Thuf 


max  of  yr  = y yjy.^  + -'jaS 


= 


(ya)2 


If 


yA  <<  1 


and 


ya  <<  1, 


it  follows  that 
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(2.1S9) 


(2.16o) 


(2.i6i) 


(2.162) 


(2.163) 


(2.16U) 


I 


max  of  yr  <<  1. 


(2.165) 


Therefore » 

= l-yr  (2.166) 

provided  (2.163)  and  fe.l6U)  are  satisfied.  This  does  not  place  any 
restrictions  on  the  relative  zone  size.  A,  with  respect  to  the  radius 
as  pointed  out  by  Tesche  [^43].  By  monentarily  letting  y = 2-n/\,  (2.163) 
and  (2.I6I4)  are  implied  by 

A « X (2.167) 

a « X.  (2.168) 

The  restriction,  (2.16?) , is  a necessity  wlicnever  a n’xnerical  teclmique 

is  forthcoming.  Accurate  reconstruction  of  the  currciit  on  an  object 
couldn't  bo  expected  if  the  zone  size  was  on  the  order  of  a wavelength. 
As  for  (2.160),  tlie  cylinder  has  ali'cady  been  assumed  thin  with  res^^ect 
to  L and  this  must  now  require  t)iat  the  cylinder  be  tliin  witli  respect 
to  wavelength.  Since  this  is  not  a steady-state  problem  but  one  involv- 
ing transients  composed  of  an  infinite  number  of  frequencies,  (2.l68) 
indicates  that  this  formulation  is  not  applicable  to  higli  frequency 
response  analysis.  In  conclusion,  the  "most"  restrictive  sentence 
encountered  is  that  this  formulation  applies  to  thin  cylinders. 
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III.  IRJMERICAL  I3iTH0DS  MID  RESIILTS 


A co;;iputer  code  has  been  vritten  to  implenent  the  equations 
developed  iti  tlie  previous  section  and  thereby  determine  the  natural 
resonances,  natural  node  vectors,  normalization  coefficients,  and 
transient  current  response.  The  cylindrical  scattcrer  is  described  by 
a general  length,  radius,  and  height  above  the  ground  plane.  Illumi- 
nation is  provided  by  a step- function  plane  wave  with  ai'bitrary  angles 
of  incidence.  Brief  nertion  of  some  of  the  numerical  methods  used  is 
warranted. 

In  order  to  better  adapt  the  equations  for  Zj,^  (2.112),  }Ip_^(2. 103) , 

and  Hi  (2.110)  to  numerical  evaluation,  their  form  can  be  altered, 
mn 

After  suppressing  the  functional  dependence  and  liighlighting  the 
matrix  notation,  these  equations  are 


_ v/l-(--/i;)-[]lr.(ia+l  .n)  - (Y2A2+2)!!o(n,n)  + HQ(n-l,n) 


‘•mn  2TiA2 


-Hj^(ra+l,n)  + (Y^A^+2)iIi(n,n)  - iri(ni-l,n)] 


(3.1) 


/-(n-l/2)A  r2^ 

Ho(m,n)  = / / 

*^(n-3/2)A  J ^ 


/h)  cos  n 


-Yri 

O '*"  I t 

d^.o  dz 


(3.2) 


/■(n-l/2)A  /•2ir 

Hi(m,n)  = / ^ 

•^(n-3/2)A  ^0  ^ ~ 


(a/h)cos  4>i 


r2 


(3.3) 


wh  »rf 
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(3.M 


ri  = [d2  sin^  ^ + ((ni-l)A-2')2]l/2 

r2  = [(ci2  + 8ah)  sin2^  + + ( (m-l)A-2' . (3.5) 

Knowlcdce  of  the  syr^netries  occuring  in  tliese  matrix  equations 
can  significantly  reduce  the  amount  of  computer  time  required 
for  this  numerical  solution.  In  (3.2),  let 

u 3 2'-(n-3/2)A  , (3.6) 


such  that 


“ r £'  r4-(- 


a/h)cos^o  r 


e d(J)o'  du 


(3.7) 


with 


r = (d^  sin2  + (A(in-n+l/2)-u)2]l/2 


(3.0) 


Obviously, 


Ho(m,n)  = Ho(r.i+i,n+i) 

i = 0,1,2,3... 


(3.9) 


since  r is  dcteiTiined  by  the  difference  between  m and  n - not  their 
actual  values.  Note  also  that  a similar  change  of  variables  in  (3.3) 
leads  to  the  same  conclusion  v?ith  regard  to  Hj^.  On  relating  this 
symmetry  to  (3.1),  it  is  apparent  that 

Zron  “ ^in+i  n+i  ^ ~ 0,1,2....  . (3.10) 

Thus,  every  element  in  any  diagonal  of  the  impedance  matrix,  Z(s),  is 
Idertical.  This  reduces  the  number  of  matrix  elements,  have 
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2 

to  be  calculated  fron  (N-2)  to  only  2N-5.  Furthermore,  upon  numerical 
evaluation  it  was  found  that 

HQ(m,n)  = HQ{n,m)  (3.11) 

and  likewise  that  this  symmetry  existed  in  Through  (3.1),  it 
follows  that 

z(m,n)  = z(n,m).  (3.12) 

Proving  (3.11)  analytically  was  pursued  at  length,  but  not  established. 
However,  this  synnetry  was  used  to  conserve  computer  time  in  calcu- 
lating z(m,n).  These  simplifications  require  that  only  HQ(m,n)  and 
H^(m,p)  for  m = 1 and  n = 2,3,...N-1  be  determined  along  with  T (2.125) 
in  order  to  form  the  entire  impedance  matrix. 

Since  both  single  and  double  numerical  integrations  were  necessary, 
a n-dimensional  Romberg  integration  routine  [47]  v;as  modified  for 
dealing  with  complex  functions. 

Other  researchers  [1-3]  have  considered  the  interactions  of  thin 
cylinders  with  an  electromagnetic  pulse  over  perfectly  conducting 
grounds.  The  "thin-wire"  assumptions  and  approximations  have  been 
incorporated  into  an  SEM  analysis.  One  of  several  computer  programs 
written  by  Shumpert  [1],  [3]  was  secured  and  data  obtained  concurrently 
with  the  computer  code  written  by  this  author.  In  this  manner,  the 
number  of  zones  used,  accuracy  requested,  and  geometrical  parameters 
were  insured  to  be  identical.  The  numerical  results  produced  by 
Shumpert's  [1]  computer  program,  shall  be  associated  with  the  title 
"approximate"  kernel,  in  the  sense  that  the  "thin-wire"  approximations 
were  incorporated;  and  in  particular,  the  current  was  approximated  by 
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filaments  of  current  on  the  cylinder  axis.  Data  found  through  equations 
developed  in  this  paper  shall  be  associated  with  the  label  "exact" 
kernel,  for  obvious  reasons. 

Before  presenting  data  related  to  the  cylindrical  scatterer  over 
the  ground  plane,  a brief  analysis  of  the  cylinder  in  free  space  is 
appropriate.  By  setting  the  image  terms  (2.110)  equal  to  zero  and 
letting  h » a,  the  impedance  matrix  (2.112)  reduces  to  one  character- 
istic of  a cylindrical  scatterer  in  free  space  with  uniform  circumferen- 
tial variation  of  the  axial  current.  The  kernel  is  still  "exact",  in  the 
sense  that  the  current  resides  on  the  cylinder  surface.  Figure  3-1 
illustrates  the  position  of  some  of  the  complex  natural  resonances, 
or  singularities,  for  both  the  "exact"  and  "approximate"  kernel.  Only  the 
second  quadrant  is  sho^sTt  - the  complex  conjugates  appearing  in  the  third 
quadrant.  As  noted  by  Tesche  [37],  the  singularities  occur  in  layers 
and  can  be  described  by  two  subscripts,  sj^j^:  denoting  the  layer 

of  the  pole  and  "n"  referring  to  the  pole  within  that  layer.  The 
natural  frequencies  are  presented  as  having  mode  distributions  that  are 
either  even  or  odd  with  respect  to  the  scatterer  center.  Figure  3-1 
illustrates  the  soundness  of  the  third  "thin-\dre"  approximation  - 
currents  might  as  well  be  approximated  by  filaments  of  current  on  the 
cylinder  axis,  since  singularity  locations  for  the  "exact"  and  "approxi- 
mate" kernels  occur  at  essentially  the  same  points.  Computer  time  required 
to  locate  the  singularities  through  the  "exact"  kernel  was  roughly  twice 
that  necessary  with  Shumpert's  [l]  program.  Table  3-1  gives  the  numerical 
data  depicted  in  Figure  3-1.  Figure  3-2  points  out  the  well-known  fact 
that  a>  the  cylinder  becomes  fatter,  the  singularities  must  move  to  the 
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left  and  down  [37].  From  circuit  theory,  this  movement  points  to  a 
decreased  quality  factor,  Q,  [48],  which  is  expected  as  the  radius  of 
the  cylinder  becomes  larger  with  respect  to  its  length.  Data  shown  in 
Figure  3-2  for  the  case  of  L/a  = 10  can  only  be  considered  as  a "first-cut" 
approximation  because  of  the  thin  cylinder  restrictions  that  have  been 
imposed  upon  the  "exact"  kernel.  Nevertheless,  one  would  expect  tlie 
results,  displayed  in  Figure  3-2,  for  the  "exact"  kernel  to  be  more 
accurate  as  the  cylinder  approaches  L/a  = 10  than  the  trajectory  produced 

by  the  "approximate"kernel  Table  3-2,  presenting  Figure  3-2  in  tabular' 
form,  shows  the  divergence  of  correlation  between  the  "approximate"  and 
"exact"  kernels  as  the  cylinder  becomes  fatter.  Note  that  in  this  figure 
and  table,  the  commonly  used  shape  parameter  is  defined  by 

f2  = 2 ln(L/a).  (3.13) 

Suppose  a cylinder  v^ith  a shape  parameter  of  10.6  (L/a  = 200)  is  now 
placed  in  proximity  to  and  parallel  to  a perfectly  conducting  ground  plane 
of  infinite  extent.  As  tlie  cylinder  approaches  the  ground  plane,  tlie 
circumferential  variation  of  the  current  becomes  nonuniform  and  one 
would  expect  the  singularities  nearest  the  imaginary  axis  to  approach 
the  theoretical  resonances  of  the  ideal  two-wire  transmission  line.  As 
displayed  in  Figure  3-3  the  singularity  associated  with  the  first  res- 
onance of  the  scattcrer  itself,  moves  toward  its  expected  destination, 

uL/c  = ir.  Both  the  "approximate"  and  the  "exact"  kernel  yield  the  same 
trajectory  until  h/L  = . 1 or  when  the  scat^rer  is  approximately  twer:^  radii 
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away  from  the  ground  plane.  The  "exact"  kernel  continues  on  its  upward 
arc  toward  the  first  resonance  of  the  transmission  line  as  the  "approxi- 
mate" kernel  falls  rapidly  toward  the  origin.  The  increasing  difference 
between  the  s^^  singularity  locations  predicted  by  the  two  different 
formulations  is  highlighted  in  Table  3.3,  With  the  cylinder  two  radii 
away  from  the  ground  plane,  the  "approximate"  kernel  was  ill  conditioned 
and  the  singularity  could  not  be  located.  Figures  3-4,  3-5,  3-6,  and  3-7 
present  similar  trajectories  of  the  s^^^  singularity  for  L/a  = JOO, 

30,  20,  and  10  respectively.  Tables  3-4,  3-5,  3-6,  and  3-7  exhibit  the 
numerical  data  observ'cd  in  the  figures  of  the  same  number.  It  should  be 
pointed  out  that  in  each  of  these  trajectories,  the  "exact"  and  "approximate" 
kernel  data  seem  to  begin  diverging  when  the  spacing  (h/L)  is  approximately 
equal  to  ten.  Tliis  fact  seems  to  be  rather  independent  of  the  particular 
L/a  ratio.  The  trajectories  of  the  singularities  associated  with  the 
second  resonance  of  the  scatterer  itself  with  L/a  = 200,  100,  30,  20,  and 


10  are  found  in  Figures  3-8  through  3-12.  As  was  the  case  with  the 
trajectories,  the  path  traversed  by  s^^^  these  figures  is  predicted  by 
both  kernels,  until  some  point  of  divergence.  Tlie  exact  kernel  natural 
resonance  continues  to  advance  towards  the  imaginary  axis  while  the 
approximate  kernel  breaks  downward.  Note  that  when  L/a  is  equal  to  200, 
the  point  of  departure  between  the  two  formulations  occurs  when  the 
cylinder  axis  is  slightly  cl^oser_JCdiafr'L7irO'Tb~the  ground  plane,  or  when 
its  axis  is  ^^-raOTi  away.  With  L/a  = 100,  the  last  point  of  agreement 
in  singularity  locations  between  the  exact  and  approximate  kernels  occurs 
when  h/L  = .1  and  with  the  cylinder  as  close  as  10  radii  away  from  the 
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ground  plane.  The  ratio  of  height  to  radius  is  apparently  not  as 
critical  as  the  ratio  of  height  to  cylinder  length.  Tables  3-8  through 
3-12  present  in  tabular  form  the  trajectories  of  Sj^2  corresponding 
Figures  3-8  through  3-12. 

The  real  and  imaginary  parts  of  the  normalized  natural  modes 
associated  with  the  first  three  resonances  of  the  scattcrer  for  L/a  = 30 
and  h/L  = 0.05  are  found  in  Figure  3-13.  It  may  be  pointed  out  that  these 
distributions  arc  essentially  indentical  to  those  found  previously  for 
scatterers  isolated  in  free  space. 

With  the  coupling  coefficient  as  defined  in  cq.  (2.146),  calculations 
of  this  coefficient  versus  the  angle  6^  were  made  for  the  case  L/a  = 30 
and  h/L  = 0.05.  These  calculations  are  shov/n  in  Figure  3-14.  (Note  the 
angle  0^  was  held  constant  at  0^  = 180°.)  As  would  be  expected,  then 
and  ni3  peak  for  broadside  incidence.  However  it  is  more  difficult 
to  interpret  the  behavior  of  n|2*  figure  3-15  presents  similar  data  for 
for  the  case  L/a  = 20  and  h/L  = 0.075. 

With  the  thin  cylindrical  scatterer  far  removed  from  the  ground 
plane,  let  the  direction  of  propagation  of  the  incident  field  be  such 
that  6^  = 30°  and  ({>^  = 180°  (see  Fig.  2-2).  The  current  on  the 
cylinder  at  z/L  = .75  is  presented  in  Figure  3-16  for  both  the  exact 
and  approximate  kernel.  Data  for  tlie  approximate  kernel  was  found  by 
Umashankar  [2],  [3]  and  is  made  note  of  only  to  illustrate  correlation 
between  the  two  kernels  when  the  circurifcrcntial  variation  of  the  axial 
current  is  uniform.  All  of  the  current  plots  in  this  note  were  constructed 
using  tlie  first  three  singularities  nearest  the  imaginary  axis,  which  is 
a valid  approach  [1-3]  when  considering  "late  time",  low  frequency 
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response.  Obviously,  the  utility  of  the  exact  kernel  formulation  is  with 
the  scattcrer  near  the  ground  plane.  With  the  scatterer  1/20  of  its 
length  away  from  the  cylinder  broadside  (0^  = 90°)  and  from  above 
(0^  = 180“).  As  shown  in  Figure  3-17,  the  transient  current  on  the 
cylinder  at  three  different  positions  is  presented  for  the  case  when 
L/a  = 30,  ^ = 0.05,  6^  = 90°,  and  (})^=  1S0°.  As  indicated,  the  low  damping 
constant  results  in  considerable  ringing  of  this  current.  A similar  tran- 
sient current  was  obtained  for  the  case  L/a  = 20,  ^/a  = 0.075,  = 90°, 

and  (J)^  = 180*^.  Again,  considerable  ringing  is  cvidneced  as  indicated  by 
the  close  prox'.  nity  of  the  poles  to  the  axis  in  the  trajectory  curves. 

Since  both  time  histories  presented  in  tlie  previous  two  figures  represented 
behavior  for  broadside  incidence  from  above  the  scatterer,  another  more 
general  case  is  included  in  Figure  19  where  L/a  = 30,  = 0.05,  o'  - 30°, 

((>^  = 180°,  and  /L  = 0.25,  0.5,  0.75. 


r 
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Due  to  the  small  difference  between  singularity  loeatioii 

for  the  exact  and  approximate  kernels  witli  respect 

to  ernph  dimensions,  they  are  shown  at  the 

same  point.  S15 j 


= 200  (2=10.6) 


EVEN  MODES 
ODD  MODES 


Table  3-1.  Location  of  the  Singularities  for  the  Free  Space  Case  where 
L/a  - 200.  (Note;  Approximate  and  Exact  Kernels  yield 
essentially  the  same  values  for  a very  thin  wire  scatterer.) 


SINGULARITY 

LOCATION 

^11 

-0.2575+j2.9093 

^12 

-0.3570+j5.9346 

^13 

-0.3903+j8.8286 

^14 

-0.3676+jl 1.5061 

^15 

-0.1920+jl5.9934 

^21 

-7.1023+j0.0031 

^22 

-8.0816+j4.7632 

^23 

-8.1417+j8.4358 
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EXACT  KERNEL  - EVEN  nOPES 


Figure  3-2.  Trajectories  of  the  Singularities  Associated  vith  the  First 
Tiiree  Resonances  of  the  Scatterer  Itself  for  the  Free  Space 
Case  as  a Function  of  L/a  - Approximate  and  Exact  Kernels. 
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Figure  3-3.  Trajectory  of  the  Singularity  Associated  with  ttie  First  Resonance 
of  the  Scatterer  Itself  os  a Function  of  h/L  or  h/a  with  L/a=200  - 
Approxii.iate  and  Cxact  Kernels 


61 


Table  3-3. 


'’r  ijectory  of  the  Singularity  Associated  with  the  First 
FesOi'irijico  of  the  Scattcrcr  Itself  p.s  a Rmcticn  of  h/L 
or  h/a  with  L/a  = 200  - Approximate  and  Exact  Kernels 


L/a  = 200  (Q  = 10.6) 


h/L 

h/a 

STI.'GULAHITl 
APPROXIMM’E  KEhiiEL 

LOCATION 

EXACT  KERNEL 

.75 

150 

-.3755 

+ 

j2.8313 

-.3769 

+ 

j2.8339 

.6 

120 

-.282li 

+ 

j2.780'4 

-.28)40 

+ 

j2.70ll 

.5 

100 

-.2212 

+ 

j2.7805 

-.2228 

+ 

j2.7806 

.k 

80 

-.1635 

+ 

j2.8001 

-.16)49 

+ 

J2.7998 

.3 

60 

-.1099 

+ 

J 2. 8335 

-.1113 

+ 

02.8380 

.2 

ijO 

-.0615 

+ 

j 2. 8987 

-.0626 

+ 

j 2. 8982 

.1 

20 

-.0217 

+ 

J2.9566 

-.0225 

+ 

j 2. 989)1 

.07 

Ih 

-.0121 

+ 

J 2. 9928 

-.0131 

+ 

J 3. 0257 

.06 

12 

-.0087 

+ 

j2.9256 

-.010)4 

+ 

J 3. 0391 

.05 

10 

-.00)42 

+ 

J 2. 5959 

-.0079 

+ 

j3.0532 

.01 

2 

-.0009 

+ 

J3.1150 
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Ficurc*  3-4 , 


I’rfijoctory  of  th-.'  f.ir.cularity  Aa'jociatcd  tho  First 

Kosoimncc  of  the-  Gcattcror  Itself  as  a ’.■'unction  of  h/I. 
or  h/a  v/ith  L/ii  ■■  ICO  - /'.t-iJro:<in:ite  and  Exact  Kornclr. 
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Table  3-4,  Trajectory  of  the  Sinc\ilarity  Associated  with  the  First 
Resonance  of  t)ie  Scattcrer  Itself  as  a Function  of  }i/L 
or  h/a  with  L/a  = 100  - Approximate  and  Exact  Kernels 


L/a  = 100  in 


h/L 

h/a 

.5 

50 

.1* 

UO 

.3 

30 

.2 

20 

.1 

10 

.07 

7 

.06 

6 

.05 

5 

.02 

2 

9.21) 


11  SlNGUL/vRlTY  LOCATIOII 


APPROXIMATE  KERIJEL 
-.2375  + j2.T063 
-.1767  + j2.73J40 
-.1202  + j2.780l4 

-.0689  + J2.8501 
-.0256  + J2.95lh 
-.0ll»6  H J2.955>4 
-.0103  + j2,8660 

-.OOhO  + j2.!4l0S 


EXACT  KERi'IKL 
-.21*11  + J2.706U 
-.1799  + j 2. 7331 

-.1232  + j2.7787 
-.0715  + o2. 8)176 
-.0275  + j2.95l)j 
-.0170  + j2.9933 
-.0138  + j3.C088 
-.0107  + J3.0256 
-.0035  + J3.0821 
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Table  3-5. 


Trajectory  of  the  Singularity  Associated  with  the  First 
Resonance  of  the  Scatterer  Itself  as  a function  of  h/L 
or  h/a  with  L/a  = 30  - Approximate  and  Exact  Kernels. 


L/a  = 30 

(fi  = 

6.80) 

h/L 

h/a 

SINGULARITY 

APPROX IKiATE  KERNEL 

LOCATION 

EXACT 

0.4 

12 

-0.1957  + j2.5422 

-0.2098 

0.3 

9 

-0.1387  + j2.6077 

-0.1505 

0.2 

6 

-0.0350  + j2.7000 

-0.0957 

0.1 

3 

-0.0362  + j2.8329 

-0.0454 

0.03 

2.4 

-0.0271  + j2.8500 

-0.0363 

0.06 

1.8 

-0.0277 

0.04 

1.2 

-0.0173  • 

KERNEL 
j2.5448 
j 2. 6066 
j2.6924 
j2.8208 
j2.8547 
j2.8936 
j2.9397 
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Figure  3-6.  Trajectory  of  the  Singularity  Associated  with  the  First  Resonance  of 
the  Scatterer  Itself  as  a function  of  h/L  or  h/a  with  L/a  = 20  - 
Approximate  and  Exact  Kernels 
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Table  3-6.  Trajestory  of  the  Singularity  Associated  with  the  First 
Resonance  of  ~the  Scatterer  Itself  as  a Function  of  h/L 
or  h/a  with  L/a  = 2Cr^-Approxiiiiate  and  Exact  Kerne  s. 


L/a  = 20 

X 

h/L^X 

\ 

(n  = 

5.99) 

h/a 

SINGULARITY 

LOCATION 

APPROX irWTE  KERNEL 

EXACT 

KERNEL 

0.4 

8 

-0.198 

+ 

j2.443 

-0.2211  + 

j2.454 

0.3 

6 

-0.143 

+ 

j2.516 

-0.1644  + 

j2.519 

0.2 

4 

-0.091 

+ 

j2.617 

-0.1085  + 

j2.611 

0.1 

2 

-0.041 

+ 

j2.763 

-0.0563  + 

j2.747 

0.08 

t.6 

1 -0.031 

+ 

j2.776 

-0.0463  + 

j2.784 

0.06 

1.2 

.. 

-0.0359  + 

j2.829 
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L/a  = 10  (fi=  1<.60  5) 


X Sii  SINGULARITY  - EXACT  KERIIEL 
B sii  SIIJGUL/vRITY  - APPROXIl-ATE  KERNEL 

h/L  = .2 
h/a  = 2 


h/L  = .75 
h/a  = 7.5 
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Ficure  3-7.  Trajectory  of  the  Sincularity  Associated  with  the  First 
ResonariLC  of  the  Scattercr  Itself  as  a Function  of  h/L 
or  h/a  with  L/a  = 10  - Approximate  .ind  Exact  Kernels 
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Table  3-7  Trajectory  of  the  Singularity  Associated  with  the  First 
Resonance  of  the  Scattercr  Itself  as  a Function  of  h/L 
or  h/a  with  L/a  = 10  - Approximate  and  Exact  Kernels 


L/a  =10  (n  = 4.605) 


h/L  h/a  sii  SIEGUIJ'.RITY  LOCATION 

APPROXifiATE  ker::el  exact  kernel 

.75  7.5  -.3695*  + J2.0327  -.14506  + jP.OBSO 

.6  6 -.2905  + J2.C955  -.3565  + J2. 11*91 

.5  • 5 -.21*11*  + J2.1I182  -.2995  + j2.1982 

.U  U -.1939  + J2.2ll*l*  -.21*52  + j2.2576 

.3  3 -.11*73  + J2.3001*  -.1923  + J2.3330 

.2  2 -.1007  + J2. 1*190  -.11*01  + J2. 1*31*5 
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Table  3"8  Trajectory  of  the  Siiif^larity  Acoociated  vith  the  Second 
Resonance  of  the  Scattercr  Itse] f as  a Function  of  h/L 
or  h/a  vith  L/a  = 200  - Approxi;r,ate  and  Exact  Kernels 


L/a  = 

200  (n  = 

10.6) 

h/L 

h/a 

APPROX  t; 

S22  r.ii:GuiJ’\Ri'rY 
•;ATE  KKRilEL 

I/)CATI0N 

EX/iCT  kepj:el 

.75 

150 

-.2919 

+ 

j5.9>iOl 

-.291J; 

+ 

j 5. 9*401) 

.6 

120 

-.3633 

+ 

j6.030!) 

-.3591 

+ 

j6.O3O0 

A 

80 

-.U589 

+ 

j5.83!i3 

-.U583 

+ 

j5.81)0l) 

.3 

60 

-.3385 

+ 

j 5. 7736 

-.3395 

+ 

j5.7783 

.2 

»40 

-.2OII1 

+ 

j’5. 801)1 

-.2022 

+ 

J 5. 8092 

.1 

20 

-.O77I4 

+ 

j 5. 9181) 

-.0767 

+ 

J 5. 9339 

.07 

lU 

-.0!*68 

+ 

j5.9'<90 

-.01)56 

+ 

15.99*45 

.06 

12 

-.0369 

+ 

j 5. 9207 

-.0361) 

+ 

J 6. 0178 

.05 

10 

-.025’t 

+ 

j5.760 

-.0279 

+ 

j6.0l)28 

.01 

2 

-.001)1 

+ 

16.1552 
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L/a  = 100  {u  = 9.21 ) 


J I 5.65 


-.5  -.4  -.3  , -.2  -.1 

c 

Figure  3-S.  Trajectory  of  the  Singularity  Associated  \nth  the  Second 
Resonance  of  the  Scatterer  Itself  as  a Function  of  li/L 
or  h/a  with  L/a  = 100  - Approximate  and  Exact  Kernels 
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Table  3-9.  Trajectory  of  the  Girsularity  Associated  with  the  Second 
Resonance  of  the  Scatterer  Itself  as  a l\inction  of  h/L 
or  h/a  with  L/a  = 100  - Approximate  and  Exact  Kernels 


L/a  = 

100  (il  = 

9.21) 

h/L 

h/a 

SI?  SINGULARITY 
APPROxiM'iTE  ker:;el 

LOCATION 

EXACT  KEPJIEL 

.5 

50 

-.5076  + J5.8812 

-.5807 

+ 

J5.8969 

.U 

Uo 

-.5263  + J5.7139 

-.5304 

+ 

J5.7262 

.3 

30 

-.3791  + J5.6660 

-.3848 

+ 

J 5. 67 37 

.2 

20 

-.2271  + j 5.71514 

-.2323 

+ 

j5.7201 

.1 

10 

-.0909  + j'5.B521 

-.0938 

+ 

j 5. 8664 

.07 

7 

-.0563  + j5.886o 

-.0583 

+ 

j 5. 9367 

.06 

6 

-.01414)*  + J5. 81*51 

-.0475 

+ 

j5.9673 

.05 

5 

-.0295  + 05.6256 

-.0375 

+ 

j 5. 9931 

.02 

2 

-.0124 

+ 

J6.0958 
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I 

I 

I 

Table  3-10.  Trajectory  of  the  Sinjjularivy  Associated  with  the 
I Second  Resonance  of  the  Scatterer  as  a Function 

' of  h/L  and  h/a  witli  L/a  = 50. 


L/a  = 30 

(SI  = 6.80) 

h/L 

1 

h/a 

Sj2  Sin^jularity 

Location 

Approximate  Kernel 

Exact  Kernel 

0.5 

15 

-0.8431+j5.4258 

-0.8888+j .5.5322 

0.4 

12 

-0.6403+j5.3128 

-0.6S57+j5.3652 

0.3 

9 

-0.4535+j.S.32S)9 

-0.4917+j5.3602 

0.2 

6 

-0.2826+j3.4324 

-0.313S+j5.4481 

0.1 

3 

-0.1275+j5.6261 

-0.1515+j5.6406 

0.08 

2.4 

-0.09S9+j5.6614 

-0.1212+j5.69S3 

0.06 

1.8 

-0.0646+j5.5414 

-0.0916+j.S.7669 

0.04 

1.2 

-0.0064+j4.9573 

-0.0619+j5.8585 
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Table  3-11.  Trajectory  of  the  Singularity  Associated 

with  the  Second  Resonance  of  the  Scattcrcr 
as  a Function  of  li/L  and  l-./a  with  L/a  = 20. 


L/a  = 20 

5.99) 

h/L 

h/a 

Sj2  Singularity 

Location 

Approximate  Kernel 

Exact  Kernel 

0.5 

10 

-0.884S+j5.1403 

-1.005+j5.296 

0.4 

8 

-0.6617+j5.0947 

-0.7459+j5.185 

0.3 

6 

-0.473S+j5.1454 

-0.5995+j5.20S 

0.2 

4 

-0.3032+j5.2737 

-0.3558+j5.3131 

0.1 

2 

-0. 1452+j5.4913 

-0.1851+j5.5265 

0.08 

1.6 

-0.1146+j5.5214 

-0.1509+j5.5936 

0.06 

1.2 

-0.0711+j5.2959 

-0.1131+j5.6SS7 

» 
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L/a  = 10  (0=i4.605) 

, S12  SINGUL/iRITY  - EXACT  KE’.UJEL 
QD  Sj^2  SINGULARITY  - APPROXII-IATE  XIERXEL 


5.5 


5.0 

,U)L 

c 

>1.5 


-2 


I 

-1.5 


-.5 


J >1.0 

0 


Ficurc  3-12.  Trajectory  of  the  Einr.ularU,/  Acsociated  with  the 

Second  Resonance  of  the  Scatterer  Itself  as  a Function  of 
h/L  or  h/a  with  L/a  = 10  - Approximate  and  Exact  Kernels 
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Table  3-12.  Trajectoi'y  of  the  Singularity  Associated  vith  the  Second 
Resonance  of  the  Scattercr  Itself  as  a IXinction  of  h/L 
or  h/a  vith  L/a  = 10  - Approximate  and  Exact  Kernels 


L/a  = 10 

(n  = 

J».605) 

h/L 

h/a 

s-jp  STi'OULARlTY 

L0CATI0!J 

APPROXII-L^viE  rCElil.’EL 

EXACT  KERi;EL 

.6 

6 

-1.1060  + 

-I.702I4  + j (4. 7277 

.5 

5 

81479  + j»-5388 

-1.1620  + jjt.8091 

A 

1* 

-.6571  + j >4. 5976 

-.8720  + .fl4.8323 

.3 

3 

-.>*915  + Jfl*.7092 

f 

-.652I4  + jl4.9ll4l4 

.2 

2 

-.3366  + '114.8886 

-.I4560  + J 5. 0709 
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L/a  = 200  (0=10.6) 


-1.0  I 

Ficure  3-l6.  Current  on  Cylinder  with  h/L  = .5,  L/a  = 200,  h/a  = 100,  0^  = 30®,  = lOO® 

at  z/L  = .75  - Exact  and  Approximate  Kernels 


in 


in  ma 


IV.  CONCLUSIONS 


The  use  of  an  assumed  circumferential  variation  for  the  axial 
cuiTents  induced  on  a thin  cylinder,  near  a perfectly  conducting  ground 
plane,  enhances  the  validity  of  analysis.  "Thin-wire"  approximations 
produce  reasonable  results  when  the  thin  cylindrical  scatterer  is  more 
than  one-tenth  of  its  length  above  the  ground  plane.  As  the  thin  cyl- 
inder is  brought  near  to  the  ground  plane,  the  assumption  of  only  axial 
variation  of  the  currents  begins  to  breakdowTi. 

With  the  thin  scatterer  in  free  space  or  far  removed  from  the 
ground  plane,  the  computer  time  consumed  in  locating  singularities 
th)ough  this  formulation  is  roughly  three  times  that  required  by 
Shumpert's  [l]  computer  code.  As  the  cylinder  approaches  the  ground 
plane,  the  computer  time  required  increases  markedly.  Nevertheless, 
when  the  scatterer  is  close,  one  cannot  use  "thin-v.'ire''  approximations  - 
the  exact  kernel  with  allowance  for  nonuniforra  circumferential  variation 
of  the  axial  current  is  necessary.  The  usefulness  of  this  foianulation, 
due  to  the  increased  complexity  and  calculation  time  required,  must  be 
evaluated  in  light  of  the  particular  scattering  problem  being  solved. 

As  noted  previously,  the  assumed  circumferential  variation  of  the 
axial  current  is  linked  to  a transmission  line  mode  (TLM)  approximation. 
The  addition  of  higher  order  modes  would  improve  the  accuracy  of  the 
data  obtained.  However,  the  improvement  in  accuracy  versus  the  increase 
in  time  and  effort  required  would  necessitate  careful  study. 
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APPENDIX  A 


EVALUATION  OF  A PARTICULAR  SINGULAR  IjJTEGRAL 


This  appendix  deals  with  the  evaluation  of  a particular  singular 
integral,  namely 


1 

1 + (a/h)cos  (j» 


d«|)  dz  , 


(Al.l) 


where 


r = [z^  + sin^ 


(A1.2) 


As  shown  by  Tesche  [U3],  in  order  to  accurately  treat  the  integration 
the  singularity  can  be  integrated  analj^tically.  Interchanging  the  order 
of  integration  and  rearranging  the  form  gives 


1 

1 +(a/h)cos  (J> 


e-^r 

dz  d<J> 


(A1.3) 


“Yr 

Expanding  e in  a Taylor  series  about  r=0, 

e“^  *=  1 - yr  + r£  yS  _ r3  + ....  (-1)"  r^  y”. . . . (Al.U) 
21  31  nl 

and  retaining  only  the  first  two  terms  (tlds  triincation  places  restric- 
tions on  the  problem,  which  are  mentioned  on  page  47)  yields 

e"^**  e 1 - yr  (Al.?) 
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Thus, 


> 

I 

! 


! 


T' 


I ^ 


dz 


W^)  cos  m-') 

2 

A ^ A 


-X”  r%7t^  I // ;•—/;•• 


d().  (AI.6) 


T = -yd 


-1 

Jn  i + 


(a/h)  cos  (p  1 + (a/h)  cos  ^ 

A 
2 


/•2it_^ 

^0 


r ^ 1 

L ' " 


•-  2 


d(i> 


(A1.7) 


with  r given  by  (A1.2).  Tesche  [1*3]  used  Dwight's  equation  (200.01)  [1*9] 
to  gi\'e 


. ^ dz  = 2 In 

-A  r 


2 In  (sin  <)/2). 


(A1.8) 


Thus, 


T « -yA 


1 + 


(a/h)  cos  ^ 


d<)  + 2 


r2.  ^ 

M 1 + 


(a/h)  cos  ({i 
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Provided  (a/h)^  <1,  which  is  satisfied  (see  Figure  2-l),  any  standard 
Integra?,  table  [50]  (for  instance  equation  509.  page  1*25  of  this 


reference)  gives 


f2r, 

Jn  1 + 


(a/h)cos  <) 


= 2Tr[l-(a/h)2]-l/2  . 


(Al.ll) 


Therefore, 


_ -2TtYA 


/V 


+(a/h)cos 


2 


sin^  <>/2) 


/2v 

r 


+(a/h)cos  ()> 


ln[sin  4>/2]  d(^  . 


(A1.12) 


The  second  term  on  the  right-hand  side  of  this  equation  is  non-singular 
nni^  easily  determined  by  machine  integration,  but  the  last  term  on  the 
right-hand  side  has  a singularity.  Let 


A » 2 


C In  (sin  6/2) 
1 +(a/h)cos  <^ 


From  Dwight  [U9],  page  11*0,  equation  (603-2), 


(A1.13) 


ln|sinx|  =»  -ln2  - cos  2x  - cos  l*x  - cos  6x  ...  (Al.lU) 

2 3 


for  sin  X / o. 


Let  X *=  ^72  for  $ between  0 and  2**,  such  that 


|sin  x|  “ |sin  ♦/2l  = sin  ^/2 


(AI.15) 


and 


ln(sin  ^/2)  = -ln2  - ^ 1/n  cos 

n=l 


(A1.16) 


This  allows  A to  he  expressed  as 

•2ir 

0 


•L  1 + 


ln2) 


cos  n(^ 


(a/h)cos  ^ 


d<J>  + 2 


/iiir 

Jlri 
1 +( 


(a/h)cos  <j> 


d(t> 


A = 


^ . I ^2it 

-(a/h)^jiyii  Jq  rr 


-•iiTln2 

TI 


cos  n<{) 


(a/h)cos  ^ 


(A1.17) 
dij.  (A1.18) 


Continuinc  with  the  last  term  on  the  right-hand  side, 
•2ir 

cos  n<^  _ 

I +^a/h  ^r>nr!  A ° 

■'0 


/2ir 

xf 


(a/h  )cos  (p 
2ir 
ir 


cos  n(t> 

- 1 l'i(a/h) 


cos  <) 


d^  + 


/2ir 


cos  n6 


(a/h  )cos  <J> 


d^ 


(A1.19) 


Consider  the  last  tena  on  the  right-hand  side  of  (A1.19).  The  trans- 
formation, ifr  = 3 yields 


cos  n(^ 

1 (a/h  )cos  ^ 


<i<) 


cos  nfi 

1 + (a/h  )cos  3 


d3 


Letting  6 = 3 + 2ir  results  in 


(A1.20) 


cos  nd 

1 +(a/h)cos  <1 


d^ 


cos  n(Q-2n) 

1 +(a/h)cos  (d-2n) 


do ; 


(A1.21) 


and  since  cos  (n(0  - 2it)]  = cos  nO, 


cos  n<> 

1 + (a/h  )co8  ^ 


cos  pO 

i + (a/h  )cos  6 


dO 


(A1.22) 


96 


Therefore  (A1.19)  can  he  expressed  as 


cos  n<^ 

1 +(a/h)cos  (J> 


d4> 


2 


cos  nij) 

1 +(a/h)cos  ^ 


d(>  . 


(A1.23) 


From  Dwight  [U9],  page  219,  equation  (858.536)  is 


cos  n^^ 

1 +(a/h)cos  ij) 


d(J) 


- TT  [[l-(a/h)gll/g-l]  " 
(a/h)«il-(a/h)2JJ-/^ 


(AI.2U) 


for  0 < a/h  <l,n  = 0,l 

Substituting  this  expression  into  (A1.23)  and  then  placing  (A1.23)  back 
into  (A1.17), 


A -l*ir  ln(2)[l-(a/h)2]“l/2 


^ (a/h)n(l-(a/h)i5p7^  / 


-hiT 


, I ln2  + Y'  l/n(h/a) 

>/l-(a/h)2  1 ^1 


>^[Viru7hT2-i] 


{A1.25) 


Simple  manipulation  shows  that 


l/n(h/a)’’[  '^l-(a/h)^-l]n  = l/n 


(h/a)2-i 


Let 

b “ ■^(h/a)2_i  _ h/a 


(AI.26) 


such  that 


A " 


-UlT 

Vi-(a/h)2 


ln2 


n=l 


>■] 


(AI.27) 


Since  a < h < ••  (see  Figure  2-1) 


and  by  L'  Hospital's  rule 


= li“hH 


(a/h) 

/l-(a/h)2 


Conclude  that  for  a<h<'»,-l<b<0.  Let  c = -b,  such  that 


0 < c < 1,  then 


/l-(a/h)2 


^ cn]  . 


(ai.28: 


From  [50],  page  U32,  a logarithmic  series  has  the  form  of 


ln(l+x)  = X - 1/2x2  + 1/3x3  _ i/l^x**  + ...  , (A1.29 


for  -1  < X < 1,  Thus, 


-ln{l+c)  = -c  + 1/2c2  _ 1/3c3  + l/hc^- 


(AI.30 


since  0 < c < 1.  With 


-ln(l+c)  = D ^ 
n=l 


cn 


(AI.28)  is  written  as 


A - —===  [ln2-ln(c+l)]  , 

, Vl-(a/h)2 


with 


c ■ h/a  1^1  - [l-(a/h)2]l/2  J , 


(A1.31 


(Al,32 


Using  this  expression  for  A,  originally  defined  by  (A1.13),  in  the 
equation  for  T (A1.12)  gives  the  final  result 


T « 


=:  + 2 

\/l  - (a/h)2  •'0 


M 1 + 


In 


IF 

\2 

l»a 

J + sin2  ^/2 

ku  ..  _ 

■ )ln  2 - Infh/a 

1 - (a/h)2 

1 L 

(a/h)cos  <j)  • 

d<!)  + 


- 1 


li 


(A1.33) 


This  fxujction  can  he  easily  evaluated  numerically  since  the  integreind  ic 
not  singular. 
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